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1 Introduction 



The DIS and total cross section data [1-3] from HERA have opened new av- 
enues in our understanding of strong interactions, and models [5-7] now exist 
5^ \ which provide a unified description of 7p interactions for photon virtualities 

ranging from = to = 30000 GeV^. The theoretical situation is never- 
theless not clear. 

Indeed, a wide range of data can be described [8-10] for > 2 GeV^ by 
the DGLAP evolution [11]. Several theoretical questions need however to be 
addressed in this context. Firstly, the evolution is leading twist, and hence one 
should remove higher-twist contributions from the data before one uses the 
DGLAP equation. Secondly, the evolution introduces extra singularities in the 
complex-j plane at j = 1. These singularities start to appear at the arbitrary 
factorisation scale Qq, and their re-summation leads to an essential singularity. 
No trace of it is however present in soft cross sections. Finally, DGLAP [11] 
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evolution should be replaced at small x by the BFKL re-summation [12]. 
The latter does not lead to an essential singularity in the complex-j plane, 
but unfortunately it does not seem to be stable against next-to-leading order 
corrections [13]. 

Given these problems, Donnachie and Landshoff have proposed to use the 
soft pomeron as a higher-twist background to be subtracted from the evolu- 
tion, while a new simple pole, the "hard pomeron" [4], would reproduce the 
DIS data. Furthermore, they have shown [5] that this new singularity evolves 
according to DGLAP, provided that one removes the j-plane singularities in- 
duced by DGLAP evolution, and keeps only their effect on the hard pomeron 
residue. Again the question arises whether such a new pole should be present 
in total cross sections and whether it is perturbative or not [14]. 

Finally, we have shown that in fact no new singularity is needed to reproduce 
the DIS data [6,7], provided that one assumes a logarithmic behaviour of 
cross sections as functions of u. Double or triple poles at j = 1 provide such a 
behaviour, and enable one to reproduce all soft and hard 7p data within the 
Regge region. 

How to bridge the gap between those models and QCD remains a challenge, 
as the description of the proton, being non-perturbative, remains at best ten- 
tative. However, LEP has now provided us with a variety of measurements of 
the 77 total cross sections, for on-shell photons, and of F2 for off-shell ones 
[15,16]. One may hope that this will be a good testing ground for perturbative 
QCD [17], and that these measurements will provide guidance for the QCD 
understanding of existing models. Hence it is important to build a unified 
description of all photon processes, and to explore where perturbative effects 
may manifest themselves. The natural framework for such a program is the 
"factorisation theorem" of the analytic S matrix, which relates 77, 7p and pp 
amplitudes. This theorem is based on t-channel unitarity, i.e. unitarity in the 
crossed channel, and in the case of simple poles one obtains the factorisation 
of the residues at each pole. For more general analytic structures, one obtains 
more complicated relations, which we shall spell out in Section 2. 

Furthermore, a relation between 77 and 7p processes may be of practical use as 
some of the measurements have big systematic uncertainties. As it is now well 
known [18], the LEP measurements are sensitive to the theoretical Monte Carlo 
used to unfold the data, leading to rather different conclusions concerning the 
energy dependence of the data. This problem is manifest in the case of total 
cross sections, where the unfolding constitutes the main uncertainty. In the 
case of HERA data, the measurement of the total cross section also seems to 
be affected by large uncertainties. Again, a joint study of both processes could 
help constrain the possible behaviours of these cross sections. 
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To decide whether new singularities can appear in 7J9 and 77 scattering, one 
must first recall why singularities are supposed to be universal in hadronic 
cross sections. The original argument [19,20], which we review in Section 2, 
used analytic continuation of amplitudes in the complex-j plane from one side 
of a 2-particle threshold to the other, and considered only the case of simple 
poles. It leads to the conclusion that these poles must be universal, and that 
their residues must factorise. 

We show in section 2 that one can generalise this formula for complex-j plane 
amplitudes, and obtain a formulation which is valid no matter what the singu- 
larity is, and which leads to consequences similar to factorisation. Indeed, the 
original conclusion that singularities had to be simple poles was obtained for 
reggeons having resonances on their trajectories. The pomeron case may be 
more complicated. For instance, it is possible that no resonances are present if 
the real part of the trajectory never reaches integer values because of non lin- 
earity. It is also possible to imagine situations with harder singularities, such 
as double or triple poles, or models involving simple poles colliding at t = 0. 
We want to stress here that our goal is not to decide theoretically between 
these possibilities by solving the unitarity equations, but only to provide a 
relation between various amplitudes. 

We show in the third section that such a formula may be extended to photon 
cross sections for any value of Q^, but that new singularities may be present 
in photon cross sections. If we assume as in [6,7] that no other singularity is 
present in DIS, stringent constraints come from the positivity requirement for 
77 total cross sections and F2. We show that it is possible to obtain a good fit 
to all photon data for Q"^ < 150 GeV^ by using either double- or triple-pole 
parametrisations. For total cross sections, no extra singularity seems to be 
needed, whereas we must introduce the singularities associated with the box 
diagram at high Q^. We conclude this study by outlining its consequences on 
the evolution of parton distributions and on the possibility of observing the 
BFKL pomeron. 
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2 t-channel unitarity in the hadronic case 



2.1 Elastic unitarity 



We start by considering tlie amplitudes for three related processes: 




a a b b b b 

We shall refer to the momenta of the incoming particles as p and g, and to the 
outgoing momentum of a as p' . We use the Mandelstam variables s = {p + qY 
and t = {p — p'Y. In the s channel, these diagrams describe the processes 
aa aa, ah — a6, and hh bb. The continuation of these amplitudes to the 
t channel describes the processes aa aa, aa bb, bb bb. 

We shall write Aab{l, t, rria, rrih) for the t-channel partial- wave elastic amplitude 
for the process a + b ^ a + b, and denote by the superscript (1) the physical- 
sheet amplitude, and by the superscript (2) its analytic continuation round a 
cc threshold branch point and back to the same value of t (see Fig. 1). 

Unitarity and analyticity of the amplitude then impose the following relation 
on the discontinuity through the threshold [21] ^ : 

with p,(t) = 

In the case of hadrons, this leads to a closed system of equations, as the 
initial and final states contribute themselves to the thresholds. For instance, 
we can consider three coupled equations for protons and pions, across the vrvr 
threshold (see Fig. 2): 



^ One can derive these relations from the unitarity of the S matrix [22] and we 
sketch such a derivation in Appendix 1. 
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Fig. 1. The amplitude and its continuation in the complex-t plane around a cc 
threshold. 
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Fig. 2. The amplitude and its continuation around the vrvr and pp branch-points. 
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From (2), we can write 

^(2) 



J-TTTT 
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We see that if A^J has a singularity at j = a{t), then A^^ cannot have that 
singularity, which must be shifted at some other value when one goes around 
the cut. For instance, near a simple pole at j = a, one has A^^^ ^ and 
^(2) ~ K ^ 

j-a+pTyK ' 

Equation (2) can be conveniently written in matrix form: 

Til) _ j.(2) ^ T^i)ji^T^2) 



with 



rp{i) 



A(') A^^ 

I tor t 

p-K pp 



1,2 



(5) 



and the threshold matrix 



Ptt 





To obtain the most constraining set of equations, one then writes the discon- 
tinuities of A^'^^ across the pp threshold (the amplitudes across this threshold 
are denoted by a superscript (3) and T^^^ is defined according to (5) with 
i = 3). This gives, using the same notations as before: 



y(2) _ 2^(3) ^ t^^)r^t^^) 



(6) 



with 




Putting Eqs. (4, 6) together, one then gets: 



2^(1) _ y(3) ^ 2.(1) ^^2^(2) _^ 2^(2) ji^j-is) 



(7) 



with R = R^ + Rp. Note that T^^) = r(i)\ 
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One can then solve the equation for T^^^ to obtain 
T(^) [1 -RD]= D 



(8) 



with D = T^^^ . This situation generahses that of Eq. (3): T*^^^ is a function of 

T^^^ and its singularities cannot come from singularities in the right-hand side 
of Eq. (8), because they are exactly matched by corresponding factors in the 
left-hand side. Note that this is why we needed to consider both thresholds, 
as otherwise some amplitudes would not be present in the r.h.s. of (7). 

Hence the amplitudes A^J, A^J and A^^^ have common singularities which 
can only come from zeroes of the determinant of the matrix in brackets in 
the left-hand side of (8): 

A = det{l -RD) = for j = z^, m = 1,...,N. (9) 



The determinant can then be written 

N 

A = II Am (10) 

m=l 



with Am = at j = Zm and A.^ 7^ at j = if /c 7^ m. A^j^ (where i and / 
stand for p or tt) then becomes 

A? = E(x^)+F. (Ill 

m=l ^ 



where Fu and au-m are functions of j finite at each Zm- 

This is the basis of the complex-j-plane factorisation of the amplitudes con- 
tained in T^^\ Indeed, we can write 

4^4^i-«0' = det(T«) = ^- (12) 



Multiplying both sides of (12) by A|, and using (11), this gives 



«pp; k + Ak 



O'-Kp; fc + Afc 

m^k 



A™ 



''Txp; m 



Ar 



+ AkFpp 



O'tttt; k + Ak 

a^rp^k + Ak 

mj^k 



A. 



''Trp; m 



Ar 



+ A^F^p 
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det{D) 



(^Ai...Afc_iAfc+i...A;v J 
In the hmit j —>■ Zk, ^ 0, and we get 



O'pp; k Oi-itit; k yOinp; k 



for A; = 1 



(13) 



(14) 



or equivalently 



hm 



Am) 



finite terms. 



(15) 



The relations (8) and (12) hold for all values of t above the Arrip threshold. 
Hence they must hold for all t as the amplitudes are analytic, and in particular 
they can be continued to negative t. The consequence (14) must thus remain 
true for the amplitude in the direct channel. 

We show in Appendix 1 that this relation remains the same in the case of an 
arbitrary number of elastic and inelastic thresholds: 



T^^'^(l-RD) = D 



(16) 



where T^^^ is the matrix containing all elastic amplitudes, and D is now an 
unknown matrix containing the effects of inelastic thresholds and the contri- 
bution of elastic amplitudes across the cuts. Again, the singularities of T^^^ 
must come from the zeroes Zm of 

A = det(I - RD) (17) 



and one obtains again (14) and (15). 



2.2 Properties 

(1) The factorisation relations are in general broken when one calculates the 
contribution of multiple exchanges of trajectories. However, as the above 
argument is totally general, after all the (s-channel) unitarising exchanges 
are taken into account, one must end up again with an amplitude fac- 
torising at each singularity, even if the latter is not a simple pole. 

(2) The matrix D is sensitive to the existence of thresholds associated with 
bound states, and does not know directly about quarks and gluons which 
do not enter the unitarity equations. Hence the zeroes Zm are not calcu- 
lable perturbatively. 
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(3) One could have a spurious cancellation of the singularity if some element 
of D has a zero at j = Zm- However, it is unlikely for this cancellation 
to occur for all t or for all processes. It is however possible to "hide" a 
singularity, e.g. a.t t = for pp and pp scattering. This might provide 
an explanation for the absence of an odderon pole in forward scattering 
data. 

(4) Each singularity factorises separately. Hence it does not make sense to 
consider globally factorising cross sections or amplitudes in the s, t repre- 
sentation, unless of course the amplitude can be reproduced by only one 
leading singularity. 

(5) The relations (15) lead to a definite prediction for the residues (or cou- 
plings) of the singularities above threshold t > 4m^. As no singularity 
occurs when t is continued to the physical region for the s channel pro- 
cesses, these relations remain true there. 

(6) We have mentioned that one always obtains a relation between the am- 
plitude and its complex conjugate. In fact, this is derived in Appendix 1 
as a consequence of the unitarity of the S matrix. The relations between 
amplitudes across a cut are really derived from the relation between the 
amplitude and its complex conjugate, and hold for whatever structure 
the cuts have. 



2.3 Specific examples 



Eq. (15) is usually not mentioned, and only its consequences for the residues 
of simple poles are considered. However, we have shown that it is true in 
general, and that leads to specific predictions for any singularity structure of 
the amplitudes Alj\j), e.g. for a given order of the zeroes Zm- We shall give 
here the formulae that correspond to simple, double or triple poles, which seem 
to be three possibilities emerging from fits to hadronic amplitudes at t = 
[23]. We shall refer to these relations as the t-Channel Unitarity relations. The 
case of cuts will not be explicitly considered here, although Eq. (15) holds also 
in this case. 



For isolated simple poles 

4 (1) \ ^ Ril; m 

^ii — 



one obtains the usual relations for the residues [19] 

i?22;^ = ^^^. (19) 
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If A\j^ has coinciding simple and double poles 

4? = -^ + 7-^, (20) 

j-z {j-zy 



one obtains the new relations 



D^,D22 = iD^2)^ 

DuS22 + SuD22 = 2Di2Sl2- (21) 



A« = + + (22) 



In the case of triple poles 
the relations become 



-^11-^22 — (-^] 



n2 



12 



FnD22 + 1^11^22 = 2Fi2/^i2, (23) 

-D11-D22 + 511^22 + 5*22-^11 = 2512-^12 + -0^2- 

It is worth pointing out that the double-pole relations are not the limit of the 
triple-pole relations for a vanishing triple-pole residue. Similarly, the simple- 
pole relations cannot be obtained from the double-pole ones. The reason for 
this is that the relations (12) relate the poles of order 2n to n+1, n being the 
maximal order of the pole. Hence the first double-pole relation is contained in 
the triple-pole ones, but not the second one, and the simple-pole relations are 
entirely separate. 



3 The photon case 



Photons can also be considered within the formalism of t-channel unitarity. 
However, it is unclear whether the elastic cross section can be measured, or 
even defined [24]. In practise, what is measured is the hadronic part of the 
cross section. Hence only the hadronic part of the photon wave function enters 
the measurement, and this part is not directly related to the S matrix. 

The way to circumvent this problem is to consider the photons as external 
state insertions on the hadronic S matrix. This means that they will enter 
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the unitarity equations only as external states, and will not contribute to the 
thresholds. 



In practise, the equations (12, 16) remain the same, provided that we write 
the threshold matrix R as 



R 




(24) 



This means that A = 1 — PpDpp only involves Dpp, hence singularities can 
now come from other elements of D, and det{D) can contain singularities not 
present in A, hence breaking the factorisation relations (12,15). Namely, we 
obtain 



pp 
aW 

^7P 



^(1) 

77 



D 



pp 



A 

D. 



IP 



PpDy, 



TP 



D. 



77- 



(25) 



Extra singularities can come from D^p or D^^. In the first case, the nature of 
the singularity is different in 7p and in 77, and the couphng of the singularity, 
which contains A, must be of non-perturbative origin. On the other hand, 
singularities in D^^ can be purely perturbative. 

In the DIS case, the off-shell photons can also be treated as external particles, 
and one recovers the above equations (25) and the possibility of extra singular- 
ities. Details and general formulae are given in Appendix 2, where one obtains 
an equation similar to (25), with a matrix D depending on the off-shellnesses 
of photons: 



where Qm stands for the two virtualities {Qf, Q|) of the initial states in the 
t channel, and Qout for the two virtualities (Ql, Ql) of the final states. 

We want to point out that the position of the possible new singularities can 
depend on Q^, and as the off-shell states do not enter unitarity equations, 
these singularities can be fixed in t. 
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4 Test of t-Channel Unitarity relations 



In order to test the previous equations, and to evaluate the need for new 
singularities, we shall use models that reproduce pp, 'jp and 77 cross sections. 
Previous studies [23] have shown that there are at least three broad classes of 
models that can reproduce all forward hadron and photon data. 

The general form of these parametrisations is given, for total cross sections of 
a on b, by the generic formula^ 

a':^ = iRat + H,b)/s (27) 



where Rab is the contribution of the highest meson trajectories (p, uj, a and /) 
and the rising term Hab stands for the pomeron. The first term is parametrised 
via Regge theory, and we allow the lower trajectories to be partially non- 
degenerate, i.e. we allow one intercept for the charge-even trajectories, and 
another one for the charge-odd ones [27]. Hence we use 

Rab = Y^,{Sr-±Y~,{~sr- (28) 



with s = 2z//(l GeV^), u = p.q. 

As for the pomeron term, we consider the following possibilities: 



Hab = Xab [S•]"^ 

Hab = sDab [log S + log Cab] , 

Hab = stab log^ ( ) + log [Cab] 

I \aabJ 



(29) 
(30) 

(31) 



These forms come from simple, double or triple poles in the Watson-Sommer- 
feld transform of the amplitude (see Eq. (59) of Appendix 1), in the limit of 
cos{'dt) large, so that the contribution from the integration contour vanishes, 
and that one can keep only the leading meson trajectories and the pomeron 
contribution. 



Using the asymptotic expansion of the Legendre polynomials Pi 



r(2/ + i) 
[r(/ + i)]22' 



(32) 



^ The real part of the amplitudes, when needed to fit the p parameter, is obtained 
from derivative dispersion relations [26]. 
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we obtain, by the residue theorem (see Eq. (59) of Appendix 1), the following 
contributions to the total cross section for simple, double, and triple poles: 



vlW(j,0) 



A(i)(j,0) 



A«(j,0) 



f V V (2a + l)r(2a + l) 
'^'"'-^W,) (r(a + l))22" ' ^^^^ 
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U - ay 







j log 









■^tot = 9\^\ log ( 1 '"7T.L^'\'r,o' (34) 
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(j - ")' 

^^m2; \ml) (r(a + l))22^ 



, ^Y. 2 / \ (2a + l)r(2a + 1) , , 



In the photon case, things are a little different. Looking first at the 7p ampli- 
tude with off-shell photons, we have 

|cos(^t)l = ^W- (36) 



In the on-shell limit 0, the Legendre polynomial of Eq. (32) becomes 

infinite, hence one must assume that the amplitude goes to zero in such a way 
that the limit is finite. One can take for instance 



with q-y{0) finite. Such a choice introduces a new scale that effectively replaces 
a/Q^ -y^^i^]^ q^l^Q^) in cos('i9j), and A'^-^^ with A'^-^\ In the 77 case, with ^^d 
the off-shellnesses of the two incoming photons, in order to keep the unitarity 
relations (12) for the amplitude A^^^ instead of A^^\ one needs to assume that 



and the scales q-yiQ"^) and q-y^P^) replace nip in Eqs. (33-35). 
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4-1 Regge region 



One can think of translating the minimum s/s of the pp case into a bound 
for ly/rrip, v / {rnpq^i^Q'^)) and z//g^((5^), and use the same bound in the three 
processes. Unfortunately, the situation is really more complicated because one 
cannot extract ^'^(Q^) from the data as the logz/ terms come from a com- 
bination of simple, double (and triple) poles at j = 1, which can always be 
re-shuffled among themselves. 

In the following, we shall use cuts on the natural Regge variables 2z/, and 
cos('(9f). We find that data are well reproduced in the region 



49 GeV^ 

^/2^ > 7 GeV. (40) 

For the 77 and the 7p total cross sections, as well as for the photon structure 
function where 0, cos('(?t) becomes infinite, and only the cut on 2z/ 

constrains the Regge region. 

Furthermore, in the case of one virtual photon, experimentalists measure the 
ep or the 67 cross sections. From these, one can extract a cross section for 7*^ 
or 7*7* scattering, provided one factors out a fiux factor. As is well known, 
the latter is univoquely defined only for on-shell particles: 

atot = hm ^T^F2- 41 



The fiux factor can then be modified arbitrarily, provided that the modifica- 
tions vanish as 0. This means, for instance, that we can always multiply 
the left-hand side of (41) by an arbitrary power of (1 — x). Hence one should 
in principle limit oneself to small values of x only. We find that we can obtain 
good fits in the region 

X < 0.3. (42) 



Note that in the case of two off-shell photons, experimentalists measure axr + 
o'TL + o'LT + o'LL, which is precisely the quantity entering the factorisation 
theorem. Hence no fiux factor is necessary here. 

Finally, all the residues are expected to be functions of Q^. These form factors 
are unknown, and are expected to contain higher twists. In order to check 
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factorisation, we do not want to be too dependent on these guesses. Hence we 
choose a modest region of 

Q2 < 150 GeV^ (43) 
where most of the 7*7* points he. 

We shaU consider in the next section possible extensions to a wider region. 
4-2 Factorising t- Channel Unitarity relations 

As explained above, the simple-pole singularities will factorise in the usual 
way. Note that there is no charge-odd singularity in the photon case, hence 
only the a/f lower trajectory will enter the relations. One then gets 

Y,,Y,^{P\ Q2) = Y,,{P')Y,,{Q'). (44) 

In the case of a soft-pomeron pole, one obtains similarly 

XppX^^{P^,Q^) = X^p{P'^)X^p{Q'^). (45) 

The case of multiple poles is given by Eqs. (21, 23), and can be made more 
transparent by using the forms (30, 31) which give factorisation-looking re- 
lations for the constants (but not for all the residues! — see Eqs. (21) and 
(23) -): 

fppf,,{P^ Q') = f,p{P')f,p{Q') (46) 
with f = D, C , t, d or c. 

4-3 Dataset 

For the total cross sections, we have used the updated COMPETE dataset 
used in [29], which is the same as that of [30] except for the inclusion of the 
latest ZEUS results on cross section [2] and for the inclusion of cosmic-ray 
data. 

For 7p scattering, we have used the full set of available data [1,2,28]. 

For the 77 measurements of F2 , we have used the data of [15,16], whenever 
these included the joint x and (and P^) dependence. We have not included 
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other data as they do not have points in the Regge region. Note that we have 
not taken the uncertainties in x into account, hence the values are really 
upper bounds in the 77 case. 



4-4 Previous parametrisations 

We have first considered the results using previous studies [6,7] of ''y^*^p and 
pp scattering. Making use of the t-Channel Unitarity relations (21) and (23), 
we have obtained reasonably good predictions for a^^ and F2. However, the 
formalism breaks down in the case of 7*7* scattering, because the form fac- 
tors that we used do not guarantee the positivity of the charge-even part 
of the cross sections. Re-fitting them enables one to get closer to the data, 
but the problem of negativity remains in some part of the physical region. 
Hence, at this point, the factorisation relations have one major consequence: 
the parametrisations of [6,7] are ruled out. 

We have also considered the hard pomeron fit of [5] where the charge-parity -|-1 
rising term contains two different simple poles: the soft and the hard pomeron. 
In this case, the soft pomeron residues factorise. The hard pomeron, if it is not 
present in pp cross sections (see however [14]), then comes in as a double pole 
in 77 cross sections, see Eq. (25), and produces a cross section proportional 
to u'^'' logu. Its residue will then depend on the value of A{j = ah), which is 
unknown. This means that factorisation does not say much about the hard 
pomeron contribution, which can always be arbitrarily re-scaled. It is possible 
to get good fits using these forms, but as they do not test factorisation, we 
shall not present these results here. 



4-5 New parametrisation: triple pole 

In the triple-pole case, the problem of negativity can be cured through the 
introduction of another functional form for the form factors. To convince our- 
selves that this is possible, we have fitted F2 in several bins to 

(Q2) = a(log u + bY + cu~^-^\ (47) 

From the values of a, h and c, and the t-Channel Unitarity relations, one can 
then predict the symmetric F2 (Q^, Q^). The result of this exercise is shown in 
Fig. 3. One clearly sees that there are two branches in the fit to HERA data: 
one with positive 6, and another one with negative h. Both have comparable 
X^, but one produces positive 77 cross sections, whereas the other one does 
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Fig. 3. Prediction from t-Channel Unitarity relations for the 7*7* cross section, 
including the box diagram of Appendix 3. 

not. Armed with this information, we found that the resulting form factors 
could be well approximated by the following forms: 




logd^p{Q'^)=di 

logc-yp(Q^) =co + ci 

With the form factors obtained from our fit, we have then checked that the 
7*7* cross section remains positive everywhere. 

4-6 New parametrisation: double pole 



In the case of a double pole. Fig. 3 shows that the situation is more difficult, 
as one cannot guarantee positivity. We have tried several possibilities, among 
which a further splitting of leading meson trajectories along the lines of [31], 
but found that positivity is still not guaranteed. 
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However, it is possible to obtain a good fit, positive everywhere, if one assumes 
a slightly modified version of the double pole [32]. 

Instead of taking an sDlogs term in Hat as in Eq. (30), one can consider 

Hab = SDab [Lab + log Cab] (49) 



with 



Lah(s) 



5Re[log(l + AabS') + log(l + Aabi-S) ')] 



7;\ S\ 



(50) 



Asymptotically, this gives the same form as a double pole. Furthermore, one 
can rewrite log(l+Aa6(s) ^) = 51og(s)+log(Aa{, + l/(s) The first term comes 
from a double pole at j = 1, whereas the Taylor expansion of the remaining 
term gives a series of simple poles. Hence Dab and Aab factorises according to 



D^'yi^P 1 Q )Dpp D^pi^P ^D^piQ ), 
A^^(P^ Q'')App = A^p{P^)A^p{Q^). 

We found good fits using the following form factors: 



(51) 



D^p log C = Ci 



A^p — Ai 




(52) 



4-7 Box diagram 



One new singularity may be present in 77 scattering: it is the box diagram, 
shown in Fig. 4, which couples directly two photons to quarks. This diagram 
must be present when the photons are far off-shell and pQCD applies. As we 
have explained above, it is not at all obvious that it is present in the case 



18 







double 


trip. 


e 


Quantity 


AT 


, 2 

X 


, 2 / /\r 


X 


X /^v 


^2 


Q01 

ozi 


789. d24 


0.962 


o /u.oy9 


l.UoU 


^2 


DO 


57.686 


0.887 


09.9d3 




a^^ 


oZ 


19.325 


O.o04 


iO.ODO 


U.4o / 




oU 


17.546 


0.585 


2i.0DU 


U. ( ly 


^pp 


yu 


iUU.oM 


i.iib 


82.849 


0.92i 


^pp 


/IQ 


55.240 


1.127 


CO nnn 
Oo.yuu 




Ppp 


67 


93.948 


1.402 


98.545 


1.471 


Ppp 


11 


16.758 


1.523 


4.662 


0.424 


Total 


1165 


1150.500 
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Table 1 

Results of fits to a generalised double pole model and to a triple pole model, using 
the form factors of Eqs. (48) and (52). N is the number of data points for each 
sub-process. 




V a V a 



Fig. 4. The box diagram contribution. 

of total cross sections, and in fact we get better fits if we include it only for 
off-shell photons. Hence an extra perturbative singularity is needed at nonzero 

We have re-calculated it and confirm the results of [33] (see Appendix 3 ^ ) . In 
the following, we shall fix the quark masses at 



^ We want to point out that one need to calculate cJiL, cttl, c"Lr and ott separately 
and sum them to obtain the off-shell cross section. A contraction of gp^y does not re- 
sum the helicities properly [34] which probably explains the discrepancies between 
[33] and [35]. 
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"rriu = md = 0.3 GeV, 

m, = 0.5GeV, (53) 
rric = 1.5 GeV, 
rrib = 4.5 GeV. 

and the quarks are included only above threshold s = 21/ — — > Am"^. 
4-8 Results 

As we want to be able to vary the minimum value of 2u, and as the fits of 
[23] neither include the generalised dipole nor use 21/ as the energy variable, 
we have re-fitted the pp and pp cross sections and p parameter together with 
those for and ^(*)^(*), and imposed factorisation of the residues. 




Vs (GeV) Vs (GeV) 

Fig. 5. Fits to the total cross sections and to the p parameters. 

We show in Table 1 the P^r data point and the number of points for each 
process. We see that one obtains a very good global fo^^ both models. It 
is well known [23] that the partial fo^^ '^pp ^"^^ Ppp never reach low val- 
ues, presumably because of the presence of contradictory data. We show the 
corresponding curves in Fig. 5. 

The values of the parameters are given in Table 2 for the triple-pole and the 
double-pole cases, and the form factors are plotted in Fig. 6. 

We see that the intercepts of the leading meson trajectories are close, in fact 
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Table 2 

Parameters (in natural units) of the global fits. 

closer than those of [23]. This is due to the smaller energy region, and to the 
much larger influence of photon data on a+. 

It may also be noted, in the double-pole case, that the parameter 5 is close 
to the hard pomeron intercept of [5]. At high Q^, because the form factor A 
falls off, the logarithm starts looking like a power of 2z/, and somehow mimics 
a simple pole. It may thus be thought of as a unitarized version of the hard 
pomeron, which would in fact apply to hard and soft scatterings. 

In the triple-pole case, this is accomplished by a different mechanism: the scale 
of the logarithm is a rapidly falling function of Q^, and hence the log term 
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Fig. 6. Form factors of the triple pole (left) and double pole (right) parametrisations. 

becomes relatively more important at high Q^. Interestingly, when one writes 
the triple-pole parametrisation as a function of x and Q^, one obtains only 
very small powers (of the order of 0.1) of Q^, which do not contain any higher 
twists, contrarily to the soft pomeron of [5]. 



4-9 Total 'fp and 'yy cross sections 



We see from Table 1 that one obtains an excellent for ^/2z/ > 7 GeV, for 
a total of 62 points. The curves are shown in Fig. 5. The fit can in fact be 
continued to \/2u = 2 GeV, with a x^/point of 0.74 for 219 points. 

We have checked that adding the box diagram leads to a slight degradation of 
the fit, whether one fits the total cross sections alone or with all other data. As 
the contribution of the box is calculated perturbatively, one might object that 
one cannot use the result down to = 0, and that only the u dependence 
should be kept. Hence we have also tried to add an extra term, proportional 
to logu/u in the total cross section, but found that the fit prefers to set the 
proportionality constant to zero. Hence it seems that this singularity is not 
needed at = = 0. However, because of large uncertainties in the data, 
it is not possible to rule it out altogether. 

Similarly, we do not find the need to introduce any new rising contribution. 
However, it is clear in view of the large uncertainties that it is not possible 
to rule out completely such a possibility. In fact, our fit prefers the 77 data 
unfolded with PHOJET [36], which rise more slowly than those unfolded with 
PYTHIA [37]. Interestingly, as we reproduce both HERA and LEP data, for 
nonzero, it is not true that an extrapolation of the nonzero data leads 
to a higher estimate of the 7p and 77 cross sections. Our fit can be considered 
as an explicit example for which such an extrapolation leads to a cross section 
on the lower side of the experimental errors. 
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Fig. 7. Fits to -Pf in the high-Q^ region. We sho"w only graphs for which there are 
more than 6 experimental points. 



The fit to F2 has quite a good as well. We have checked that one can easily 
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10^ 10' 10'^ 10^ 10"* 10''' lO'" 10"* 10'^ lO"" 10' 10" 




Fig. 8. Fits to F2 in the low-Q^ region. We show only graphs for which there are 
more than 6 experimental points, as well as the lowest-Q^ ones. The curves and the 
data are as in Fig. 7. 

extend it to Q'^ ^ 400 GeV^ for the triple pole, and to ^ 800 GeV^ in 
the double-pole case. It is interesting that one cannot go as high as in ref. [7]. 
This can be attributed either to too simple a choice for the form factors, or to 
the onset of perturbative evolution. 
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Figs. 7 and 8 show the fit for the most populated bins. As pointed out 
before, we see that our fits do reproduce the low-Q^ region quite well, but 
predict total cross sections on the lower side of the error bands. Hence the 
extrapolation to = of DIS data does not require a hard pomeron. 

J^.n Fits to F1 
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Fig. 9. Fits to F^. The curves are as in Fig. 7. The data are from [15,16]. 

As the number of data points is dominated by and 7p data, the fit to 77 data 
is really a test of the t-Channel Unitarity relations. As we explained above, 
the strongest constraint comes from the positivity of the 7*7* cross section, 
which is not guaranteed by the t-Channel Unitarity relations in the case of 
multiple poles. As Fig. 9, 10 and 11 show, one obtains a good description of 
the points within the Regge region. 

Here, we have observed that the quality of the fit improves if we add the box 
diagram for nonzero and . There is no need however to include other 
singularities, such as a hard pomeron or a perturbative one. 



25 



S.SGeV^ 









10"^ 10"^ lO"' 10"^ 10'^ lO"' 



X X 

Fig. 10. Fits to F2 for = Q^. The curves are as in Fig. 7. The data are from [15]. 



p2 = 3.7 GeV^ 
Q^ = 120 GeV^ 
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X 



Fig. 11. Fits to F2 for nonzero asymmetric values of and . The curves are as 
in Fig. 5. The data are from [15]. 

For 7^ and = 0, the box diagram makes Uttle difference in the double- 
pole case, but does reduce the appreciably in the triple-pole case. We have 
included it in both cases. 



5 Conclusion 



We have shown in this paper that t-channel unitarity can be used to map the 
regions where new singularities, be they of perturbative or non-perturbative 
origin, can occur. Indeed, we have seen that although hadronic singularities 
must be universal, this is certainly not the case in a^p^hadrons, F2 or F2 , as the 
photons enter only as external particles through an insertion in the hadronic 
cross section. 

We have shown however that up to^ = 150 GeV^, the data do not call for 

^ The region we have considered excludes the highest-Q^ 77* points from OPAL. 
For the point which falls in the Regge region, at = 0, = 780 GeV^ and 
X = 0.275, the experimental value is 0.93 ± 0.16, the extrapolation of the double- 
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the existence of new singularities, except perhaps the box diagram. 



For off-shell photons, our fits are rather surprising as the standard claim is that 
the perturbative evolution sets in quite early. This evolution is indeed allowed 
by t-channel unitarity constraints: it is possible to have extra singularities in 
off-shell photon cross sections, which are built on top of the non-perturbative 
singularities. But it seems that Regge parametrisations can be extended quite 
high in without the need for these new singularities. 

Finally, as the BFKL singularity is purely perturbative (the position of the 
singularity and the form factor come from pQCD), it can manifest itself only 
in 7*7*, but we have seen that there is no definite need for such a singularity 
in present data. 



Appendix 1: general proof of t-channel unitarity relations 

Assuming that m;, is the lowest hadronic mass, we know that aa — *• aa, ah —>■ ah 
and hh — > hh have thresholds for t > 4m^ > 4m^. In general, if t is large enough, 
there are many possible intermediate states for each process under considera- 
tion, which we must take into account to write the unitarity relations. These 
states can be grouped into subsets which have the same quantum numbers, 
and for which one can derive factorisation. 

Starting with the unitarity of the S matrix: 



(54) 



and setting S* = E + iSc, we obtain 



(55) 



One can define the invariant amplitude Tif by the matrix elements 



<f\S,\z>= {2iTyS\pf-p,)T,j. 



(56) 



Eq. (55) then becomes the following at the amplitude level: 



T,f-T^f = Cs{T,T^). 



(57) 



pole fit predicts 0.71, while the triple-pole prediction is 0.74. 
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We define the Cg operator as the following convolution: 

Cs{T\ T) = CsiT, Tt) = 2z 5^ 1 dPS T.^tI^ (5J 



where k labels the possible intermediate on-shell n-particle states in the t 
channel, which can differ by the number and nature of produced particles, 
and dPS represents the differential n-particle Lorentz-invariant phase space 
associated with these states. 

If the particles are massive, we can enumerate these open channels and assume 
that k runs from 1 to + 2, the number of open channel depending on the 
value of t, i.e. of the energy in that channel. In particular, we shall find in 
this set of states the aa and hb intermediate states to which we respectively 
assign the labels k = 1, 2. Note that in general the label k does not refer to 
the number of particles in the intermediate state, and that k can stand for 
particles different from a and b. So in general the amplitude Tkm represents 
the following process: 

state k 




We can now imagine that we split the amplitude into charge-parity -|-1 and 
charge-parity —1 parts, and perform a Watson-Sommerfeld transform 

r 2,1 -\- \ 

r±(z.,t) = 167r2yt^/PK-cos(^0)^-^Tat(/,i) (lie-^'^'). (59) 
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with V = p.q. (In the following, we shall only consider the charge-parity +1 
part of the amplitudes without carrying the superscript +.) After continu- 
ing this relation to complex / = j, we deform the contour of integration so 
that only the singularities of T(j, t) will contribute. All amplitudes become 
functions of j, and the operator Cs changes to C, which has the following 
properties: 

• It is associative and distributive 

C{aAi + f3A2, B) = aC{Ai, B) + (3C{A2, B). (60) 

• In the case of 2-particle intermediate states /c, the form of C is particularly 
simple: 

C2{T\ T) = pkT..kTlf = {TRT^),f (61) 



/ ^ 4:771 

with pk = 2i^—^, and Rkm = pAm- 

To proceed further, we shall represent the T matrix in the following form, for 
k<N + 2: 



To (2 X 2) 


Tu{2 X N) ^ 


Ti{Nx2) 


Tr{N X N) ^ 



(62) 



where we have indicated the dimensions of the sub-matrices in parenthesis. Tq 
contains the elastic amplitudes (^j-/\ /=1, 2), the upper matrix contains 
the inelastic amplitudes i = l,2— >/c>2, and the lower matrix Ti the inelastic 
amplitudes k > 2 ^ i = 1^ 2. stands for the rest of the amplitudes k ^ m, 
with k and m > 2. 

The system (57) can then be written: 



To-Tj = Toi?Tj + C(r„,Tt), 
T^-T^ = ToRTI + CiT^,Tj), 
Ti-Tl = TiRTl + C{Tr,Tl), 

Tr-T} = C{TuTi) + C{TrX)- 



(63) 
(64) 
(65) 
(66) 



To derive factorisation, it is enough to consider the first two relations (63, 
64). We assume that the second equation can be solved by a series expansion, 
yielding 

T„ = M + TqRM (67) 
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with M the solution of M = T/ + C(M, T^): 



M = T} + C{Tl Tl) + C{C{Tl T}),Tl) + ... 



(6J 



We can put this form into Eq. (63), which then gives 



To{l-RD) = D 



(69) 



with 



D 



t^ + c{mX] 



(70) 



Then we can repeat the argument given in the main text after Eq. (8), leading 
to the factorisation relation (14) near each singularity. 



Appendix 2: unitarity constraints for off-shell photons 



The virtual photons must not be included in the intermediate states of 
Eq. (57). For hadronic final states, and in the one-photon approximation, 
the photons can be thought of as external particles which get inserted into the 
hadronic cross section. 



In this case, we want to indicate explicitly whether the external legs of the 
2 ^ 2, 2 ^ 77, and n ^ 2 amplitudes are off-shell or not. We introduce the 
notations To{Qin,Qout), Tu{Qin) and Ti{Qout), where Qm stands for the two 
virtualities (Qi, Q2) of the initial states in the t channel, and Qout for the two 
virtualities {Ql,Ql) of the final states, and we write Qm = or Qout = in 
the case of on-shell states, and the relations (57) can be visualised as follows: 
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The system of equations (63-66) then becomes: 

To{Qin, Qout) ~ 2^0 (Qin, Qout) = T{){Qin, 0)-RTq (0, Qout) 

+ C{TM,^)XiQout)), (71) 

Tu{Q^n) - T/(g„) =To(Q,„, 0)i?T/(0) + C(T„(Q.„), T,t), (72) 

TiiQout) - TliQout) =Ti{Qout)RT^{0, 0) + C(T„ Tj), (73) 

Tr-T} = C{Ti{0),T^{0))+C{Tr,T}). (74) 

The resolution of the system proceeds as in Appendix 1 with the elimination of 

Tui^Q in) ■ 

Tu{Q^n) = M{Q,n) + T^{Qin. 0)i?M(0) (75) 

with M the solution of M{Q) = T^{Q) + C{M{Q), T}): 

M{Q.n) = Ti\Q,n) + C(T/(g,„),T,t) + C{T^{Qi^),Tj),Tj) + ... (76) 

The first equation however now gives 

ToiQin, Qout) = D{Qin, Qout) + T^^Qin, 0)RD{0, Qout) (77) 

with 

D{Qin, Qout) = T^iQin, Qout) + C (M (Q.n) , T^^Q out)) ■ (78) 

For DIS, we consider Qout = and Ql = Q^ = Q^ = —q^. (Note that the same 
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kind of relations and conclusions would hold for off-forward parton distribution 
functions). This gives us 

ToiQin, 0)(I - RD{0, 0)) = D{Qin, 0). (79) 

Hence we see that all the on-shell singularities must be present in the off-shell 
case, due to the factor (E — RD{0,0)), but we can have new ones coming 
from the singularities of D(Qj„,0). These singularities can be of perturbative 
origin {e.g. the singularities generated by the DGLAP evolution) but their 
coupling will depend on the threshold matrix R, and hence they must know 
about hadronic masses, or in other words they are not directly accessible by 
perturbation theory. 

In the case of 7*7* scattering, we take Qm = and Qout = = — p^, and 
Eq. (77) gives 

ToiQin, Qout) = D{Qin, Qout) H ^j^jjf^ i^^) 

This shows that the DIS singularities will again be present, either through A, 
or through extra singularities present in DIS (in which case their order will be 
different in 77 scattering, at least for Qj„ = Qout)- 

It is also possible to have extra singularities purely from D{Qin, Qout)- A priori 
these could be independent from the threshold matrix R, and hence be of 
purely perturbative origin [e.g. 7*7* — ^ qq or the BFKL pomeron coupled to 
photons through a perturbative impact factor). 



Appendix 3: the box diagram 



We have re-calculated the contribution of the box diagram of Fig. (4.7), and 
confirm the results of [33]. Our results can be recast in the following form, 
which may be more transparent in the present context: 

We use xi = P^/ {2iy) and X2 = Q^/ (2i^), with z/ = p.q, which give 



1 - Xi- X2 



Q-' = y^^^ (82) 

I - Xi- X2 

with w"^ = s. We set 



2 
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li- 



ra 



2m?{\ — Xi — X2) 



r = 1 — 4xiX2, 
6 = — Xi — X2 + 1, 
Sf, = 6 - 2fi. 

The cross sections then take the form 



6t{26xiX2 + r/i)r2 r 



^ log(p) 



(83) 

(84) 
(85) 
(86) 



^7) 



which gives 



P 



55^T + 5f,T 



The cross sections then are built from: 



TiTT = 4(5xiX2 [2x1X2 (x^ + X2 — 1 + 2xi + 2x2) 

- 12x1X2 - 2(xi + X2) + 2(x2 + xi) - 1] 

- r/i(2xi - 1)^(2x2 - 1)^ - 25/iV^ 
Att = '26fiT - 2/iV2 

+ [8x^X2 (x^ + X2) + 16x'^X2 — 16x^X2 (Xi + X2) 

-4xiX2(xl^ + X2) + 16x1X2 + -2(xi + X2) + 2(xi + X2) + 1], 
Stl = /ir5x2[(6xi + 1 + 2x2X1 - 6x1) 

+ 25xi((2xi + l)x2 + (2x2 + - 6x1X2)], 
Atl = — X2[2(5xi(2x2 — 1 — 2x1X2 — 2x1X2 — 2x1X2 + 2xi) 

+ /ir(2xi + 1 - 2x2X1 + xi)], 
Tjlt = ^tl{xi ^ X2), 
Alt = Atl{xi ^ X2), 
Sll = -25^x1X2(35x1X2 + /ir), 
All = -5^X1X2(2x1X2 + 1). 
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